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Measurements have been conducted at the University of Virginia Supersonic Combustion 
Facility of the flow in a constant area duct downstream of a Mach 2 nozzle, where the 
airflow has first been heated to approximately 1200 K. Dual-pump CARS was used to 
acquire rotational and vibrational temperatures of N 2 and O 2 at two planes in the duct at 
different downstream distances from the nozzle exit. Wall static pressures in the nozzle are 
also reported. With a flow of clean air, the vibrational temperature of N 2 freezes at close to 
the heater stagnation temperature, while the O 2 vibrational temperature is about 1000 K. 
The results are well predicted by computational fluid mechanics models employing separate 
‘‘lumped” vibrational and translational/rotational temperatures. Experimental results are 
also reported for a few percent steam addition to the air and the effect of the steam is to 
bring the flow to thermal equilibrium. 


Introduction 

Hypersonic flow simulation facilities provide high-speed, high-enthalpy flow conditions to test aerospace 
vehicles, thrusters and engines. The test gas is typically heated to avoid freestream condensation, to allow 
combustion tests in engines and/or to simulate high enthalpy “real-gas” effects such as dissociation. The test gas is 
also accelerated through converging-diverging nozzles to provide the required flow velocity. However, heating and 
expanding the test gas imparts non-ideal characteristics to the gas. For example, combustion-heated facilities can 
have lO’s of percent of water in the freestream while shock- or arc-heated facilities have a few percent of nitric 
oxide (NO) present in the test gas. This contamination is known as free stream vitiation. Other imperfect attributes 
of gases produced in hypersonic facilities include the existence of atomic oxygen and nitrogen from the dissociation 
of O 2 and N 2 respectively, or copper or iron eroding from facility electrodes or facility walls. Yet another imperfect 
characteristic that the gas can acquire in hypersonic test facilities is vibrational non-equilibrium. This phenomenon 
often occurs when a high-temperature stagnant gas is expanded suddenly through a converging-diverging nozzle to 
supersonic or hypersonic velocity. At first, the expansion occurs slowly enough (high enough collision frequency) to 
keep the rotational and translational modes of diatomic molecules in equilibrium with each other. However, as the 
gas passes down the nozzle, the expansion can become too great (too low collision frequency) and the flow speed 
too fast for vibrational modes to stay in equilibrium with the rotational and translational modes of the gas.^ The 
number of collisions required for the rotational mode to be brought into equilibrium with the translational mode is 
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much less than required to bring the vibrational mode into equilibrium, so the rotational mode is normally assumed 
to be in equilibrium. Consequently the Boltzmann distribution of population among rotational states, resulting in a 
rotational temperature, is out of equilibrium with the distribution among vibrational states, Ty. Further 
complicating this situation, different molecules (N 2 , O 2 , NO, etc.) can develop different vibrational temperatures. 
These non-equilibrium effects may be important because they can affect ignition or reaction rates in supersonic 
combustion engines tested in ground-based facilities. Ignition could be affected because molecules in excited 
vibrational states may be more reactive than those in ground vibrational states. Vibrational non-equilibrium also 
changes the static pressure measured at the exit of a nozzle because it changes the effective specific heat capacity of 
the expansion process. In order to compare computational fluid dynamics (CFD) calculations with experiments, it is 
critically important to model or to determine the magnitude of these non-ideal characteristics. Non-equilibrium 
effects can be computed but measurements should be made, if possible, to validate the computations. 

Vibrational non-equilibrium in hypersonic facilities has been measured previously in arc jets using spontaneous 
Raman spectroscopy,^ and in shock tunnels using nitric oxide planar laser induced fluorescence (NO PLIF)^ and 
broadband N 2 coherent anti-Stokes Raman spectroscopy (CARS)."^ It has been included in the calculation of facility 
freestreams and shock-layer flows.^ This phenomenon has been studied analytically or through correlations by 
numerous researchers.^’^ Vibrational non-equilibrium has been predicted to occur in electrically heated blowdown 
wind tunnels.^’^ In this paper, the dual-pump coherent anti-Stokes Raman spectroscopy (DP-CARS) technique is 
used to study the vibration-rotation non-equilibrium in the supersonic nozzle expansion at the University of Virginia 
Supersonic Combustion Facility (UVaSCF). DP-CARS is a non-intrusive, laser-based technique, originally 
developed by Lucht,^ that has been implemented by several researchers to measure temperature and species 
concentration in scramjet combustors. At the NASA Langley Research Center, the DP-CARS technique was 
successfully used to measure temperature, N 2 , O 2 and H 2 concentration in scramjet engines. 

In a previous paper, we reported DP-CARS measurements in the UVa dual-mode scramjet “Configuration 
We observed vibrational non-equilibrium in the spectra of air in the freestream of this experiment, particularly for 
the N 2 . We noted that, according to Landau-Teller theory,^ the relaxation time for pure N 2 is roughly 28 ms at the 
nozzle exit condition while the flow travels the length of the scramjet in only about 1 ms. Thus, the potential 
observation of vibrationally frozen flow was deemed plausible. 

In the present paper, we report a more detailed analysis of vibrational non-equilibrium in the UVa facility. The 
experimental data in this paper were acquired in two separate test entries. The first is the aforementioned entry in 
which measurements were made in Configuration A. The second occurred in the spring of 2012 using Configuration 
C. The facility operating conditions were the same in both configurations but the model geometry was different. 
Configuration A consisted of two sections of duct downstream of the nozzle: a Combustor section followed by an 
Extender section. For this configuration, the CARS data presented herein were obtained just downstream of the 
nozzle exit. Configuration C consisted of Isolator, Combustor, TDLAT (enables tunable diode laser absorption 
tomography measurements), and Extender sections downstream of the nozzle. For this configuration, the 
measurement plane was just downstream of the exit of the constant area Isolator section. Consequently, 
measurements at the same operating conditions allow the vibrational relaxation in the isolator duct to be studied. 
Also, for the Configuration C cases, various levels of steam were added. Addition of steam was expected to bring 
the gas flow back towards equilibrium.^^ Computations of the flowfield for both equilibrium and thermal non- 
equilibrium have been performed and are compared with the data. Our contribution is to demonstrate that vibrational 
non-equilibrium can be found in supersonic nozzles as well as hypersonic nozzles, and show that published 
modeling methods can predict it. 


Experimental Facility 

The UVaSCF is a direct-connect electrical-resistance-heated facility capable of providing continuous, clean air 
flow with a total temperature of 1200 K. In the current experiments, the facility was operated at a total pressure of 
296 kPa with an air flow rate of 0.18 kg/s issuing through a Mach 2 facility nozzle. These conditions simulate a 
scramjet vehicle traveling at Mach 5 with engine inlet compression slowing the flow to Mach 2 at the entrance to the 
engine isolator. Details of the general tunnel operation are described in Refs. 16 and 17. The experiments were 
performed with no fuel flow and no combustion in the scramjet flowpath. 

The facility has the capability of simulating the operation of a combustion-heated facility through the addition of 
steam, carbon dioxide, and supplemental oxygen that is injected into the air flow upstream of the electric heater. 
The flow rates of each gas are metered via choked orifices. Stagnation pressure and temperature are measured 
upstream of each orifice and, together with a measured orifice diameter and assumed discharge coefficient, this 
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provides sufficient information to determine the mass flow rate of the added gases. Injection upstream of the facility 
heater promotes full mixing and heating of the added gases. Only steam has been added for the experiments 
described here. 

The water vapor is supplied by an electric boiler that can generate steam continuously at pressures up to 
approximately 3500 kPa. The generated steam travels through plumbing to a three-way valve upstream of the tunnel 
that allows steady-state operation of the steam system and heater to be achieved before the steam is introduced into 
the main air flow. This is important as the air temperature at the point of steam injection must be sufficiently high to 
prevent water condensation. A combination of heat from the facility heater core and external steam lines can raise 
the incoming air temperature to 360 K to allow up to 15% H 2 O by mole fraction without condensation. 

Measurements were performed at three different levels of water vapor addition: 1.9%, 3.7%, and 5.7%. Sources 
of uncertainty in the steam delivery system include the discharge coefficient and throat area of the choked orifice as 
well as instrumentation uncertainty associated with measuring pressure and temperature. In addition, the steam 
generator itself may not maintain a perfectly constant steam pressure and temperature with time. The steam injection 
system has been previously verified and the uncertainty quantified as 3.4% of the steam flow rate.^^ 

The Mach 2 facility nozzle is a two-dimensional, symmetrical, contoured nozzle with an area ratio of 1.786. The 
contour in the vicinity of the throat consists of two circular arcs that join tangentially at the area minimum. The radii 
of the arcs downstream and upstream of the minimum, normalized by the half height of the nozzle throat (7.1 1 mm), 
are 5.38 and 9.38, respectively. Downstream of the throat, the nozzle is contoured to straighten the flow. The side 
walls of the nozzle are parallel and the exit is Ah in height and 6h between the side walls, where h = 635 mm is the 
reference length in this experiment (equal to the height of the fuel injector). 

The dual-mode scramjet consists of a length of constant area duct matched to the nozzle exit, followed by a 
diverging duct with a fuel injector. CARS measurements may be made at a series of cross stream planes using 
specially designed windows in the two sidewalls. Each incorporates a 3.18 mm wide slot that spans the entire height 
of the duct, with an airtight enclosure and offset (to avoid laser beam damage) windows at Brewster’s angle. This 
enables CARS measurements over the entire cross plane at each axial station. The measurements reported herein all 
pertain to the most upstream plane, which was located at the beginning of the combustor section just prior to the 
start of the divergence in the duct. For the purposes of this paper, the only significant difference between 
Configurations A and C is that C has an extra constant area length of duct (the “isolator”) located between the 
combustor and the nozzle, whereas A does not. Consequently, for A the CARS measurement plane is located at x/h 
= 13.8 from the nozzle throat, whereas for C it is at x/h = 55.6. 

Dual-Pump CARS Technique 

Coherent anti-Stokes Raman spectroscopy (CARS) is a non-linear spectroscopic technique in which three laser 
beams-two pump beams (green and yellow) and a Stokes beam (red)-are focused and crossed at their focal point, 
and a fourth signal laser beam (blue) is generated at the intersection through a four-wave mixing process, as 
illustrated in Fig. l(a,b). Resonances associated with Raman-active molecular transitions between rotational- 
vibrational states strengthen the CARS signal in a manner dependent on gas composition and temperature. The 
wavelength of excitation-or “Raman shift” -corresponds to the difference in energy between the pump beam photon 
and the Stokes beam (longer wavelength) photon. In the broadband CARS technique, the Stokes laser beam is 
spectrally broad whereas the pump beams are spectrally narrow, and this enables multiple resonances to be excited 
simultaneously. The signal beam is spectrally broad and carries with it the spectral signature of the gases present 
(over some range of wavelengths). In the dual-pump technique, two CARS processes (in which the roles of the 
pump beams, copj and C 0 p 2 , are exchanged) take place simultaneously, and the signals from these two processes are 
superimposed. This signal beam is dispersed by a spectrometer and spectra are compared to analytical spectral 
models to obtain information on composition and temperature. 
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Figure 1. (a) Laser beam arrangement of DP-CARS, (b) energy level diagram of CARS process 

The CARS system is the same as described in Refs. 14 and 19. An injection seeded Nd:YAG laser, frequency 
doubled to 532 nm; a narrow band dye laser (Spectra Physics PDL-2) operating at 551 nm; and a home-built 
broadband dye laser centered at 603 nm with FWHM of 10 nm, are used to generate the CARS signal. The pulse 
width of the lasers is about 8 ns and the repetition rate is 20 Hz. The green output of the Nd:YAG laser, 
approximately 900 mJ per pulse, is divided three ways: 170 mJ to pump the broadband dye laser, 415 mJ for the 
narrowband dye laser, and 280 mJ for the CARS green beam. The two dye lasers are combined using a dichroic 
mirror. The three lasers are located on a mobile cart, and the beams are relayed to the measurement volume through 
several relay mirrors. The beams are focused by a 75 cm focal length spherical lens for the superimposed dye laser 
beams and a 60 cm lens for the green beam. At the measurement volume, the three beams interact and generate the 
signal beam containing a broad CARS spectrum centered at 490 nm. The laser energies at the measurement volume 
are 63, 32 and 19 mJ for the Nd:YAG, the narrowband and the broadband beams respectively. The lens focusing the 
green beam is tilted with respect to the incoming beam so that, by introducing astigmatism, it shapes the focal spot 
at the beam crossing to a roughly 4:1 axis ratio ellipse. This helps in reducing beam steering effects in flows with 
large, unsteady density gradients.^^ A series of dichroic mirrors and filters is used to separate the signal from the 
Nd:YAG beam. The signal is then focused by a spherical lens on the input slit of a one meter spectrometer with a 
2400 grooves/mm grating. The CARS signal is dispersed by the spectrometer and recorded on a cooled, back- 
illuminated CCD camera. In the implementation at UVa, the laser beams are generated on a mobile cart that is 
located in a lab adjacent to the scramjet. Beams are relayed to the measurement plane in the facility via a series of 
dichroic mirrors and a translation system that enables the measurement point to be moved in the flow, and the signal 
is relayed by additional mirrors from the measurement volume to the detection system located on an optical bench in 
the scramjet lab.^"^ 

The collected spectra are analyzed by comparison with the Sandia CARSFT code^^ and using the fitting 
algorithm developed by Cutler.^^ The CARSFT code contains all the theoretical models for calculating the spectra, 
including the models for computing the CARS theoretical susceptibility and convolutions with laser line shape and 
instrument function. CARSFT is capable of either comparing with or fitting to an experimental spectrum, but the 
algorithm converges poorly (or not at all) with more than a couple of fitted parameters at a time. When handling 
large data sets, especially when fitting for both temperature and several species mole fractions, the Cutler code is 
used. This code interpolates from previously-generated libraries of spectra computed using CARSFT, and fits for 
experimental parameters that may vary from spectrum to spectrum (amplitude scaling and horizontal shift) as well 
as the temperature and species mole fractions. The library contains spectra at various temperatures and species 
compositions and is sparsely structured so that its size scales to a power less than the number of species plus one.^^ 
Generally, it is assumed that not all the species present in the gas being probed appear in the CARS spectrum, and 
that the nonresonant susceptibility of the species not present (the buffer gas) is constant and known. In this work, 
both N 2 and O 2 Q-branches appear in the spectrum so both major species in air are present. We use the Voigt line 
shape model^^ and take the pressure (which affects the spectra through the line shape) to be the pressure measured at 
the edge of the flow by wall static pressure taps. In the current paper, air is treated as a single species variable in the 
library (consisting of 0.79 mole fraction N 2 and 0.21 mole fraction O 2 ) while H 2 O is the buffer gas and the library 
has two variables (temperature and mole fraction of air). If only air is present, the mole fraction fits to 1.0 give or 
take an error associated with errors in the spectral modeling or assumed pressure. 

The CARSFT code was modified for this work to permit the calculation of spectra with different rotational and 
vibrational temperatures of N 2 and O 2 and comparison to experiment, but was not modified to fit for vibrational 
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temperatures. The Cutler et al.^^ fitting code was not able to fit for separate vibrational temperatures either, although 
it could employ libraries in which the vibrational temperature was set to a specified constant or set equal to the 
rotational. Specific fitting strategies are discussed in the Results section. 

Computations 



Figure 2. Computational grid 

The CFD simulations of the flowfield were performed using the LAURA solver^"^ developed at the NASA 
Langley Research Center. The conditions considered involved dry air (mass fractions N 2 and O 2 of 0.767 and 0.233 
respectively) at a stagnation pressure (pheater) and temperature {Theater) of 296 kPa and 1200 K, respectively. Two 
simulations were performed, one which assumed thermal equilibrium and one in which thermal non-equilibrium 
processes were modeled. Simulations for conditions that contained steam were not attempted, due to concerns about 
whether the modeling assumptions employed for the thermal non-equilibrium processes were valid for steam. 
Thermal non-equilibrium was modeled using a two-temperature model that assumes that translational and rotational 
energy modes (with temperatures Ttrans and Ty) for each species are in equilibrium, and that the vibrational and 
electronic energy modes for each species are in equilibrium at temperature The simulations were advanced in 
pseudo-time using a point-implicit algorithm with a Courant-Friedrichs-Lewy (CFL) number of 5. The inviscid flux 
components were evaluated using the upwind scheme of Roe,^^ and viscous flux components were evaluated with a 
centered scheme. The Menter SST model^^ was chosen to account for flow turbulence. The nozzle flowfield was 
simulated as a two-dimensional geometry, using the multi-block structured grid shown in Fig. 2. The boundary 
conditions utilized are also labeled in Fig. 2. The grid topology and resolution was based on the three-dimensional 
grid employed by researchers at North Carolina State University in their computational studies of this dual-mode 
scramjet flowpath. The surface temperature was set to 600 K along the isothermal walls to match that used in these 
previous computational efforts. 

The coordinate system used in this work is also shown in Fig. 2; x is referenced to the throat of the nozzle; y is 
referenced to the center of the duct; z is according to a right-handed system and is referenced to the center of the 
duct. 


Results 


A. Existence of Non-equilibrium 

Early analysis of Configuration A CARS spectra acquired in the freestream air of the combustor flow path, 
conducted while the Configuration A test campaign was still underway, led us to suspect “anomalies”: while the 
rotational band structure of the ground vibrational state of the N 2 and O 2 Q-branch spectra was consistent with the 
expected temperature, the first excited vibrational state intensity was too high. It was possible that the observed 
vibrational non-equilibrium could have been generated by the CARS process itself This well-known phenomenon 
is called stimulated Raman pumping (SRP). The green and red laser beams transfer ground- state population to the 
first excited vibrational state of N 2 if the laser energies are sufficiently high via a two-photon process (the left-most 
up and down arrows in the energy level diagram shown in Fig. Ib).^^’^^ To determine if SRP was affecting the 
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measurements, the energy of the green laser {copi) was varied by a factor of 34 during one run of the experiment. 
The averages of 200 spectra in each case are shown in Fig. 3 where all the spectra are normalized to equal area 
under the curve. The horizontal axis is Raman shift (cOpi - co^ in wavenumbers. If SRP were the cause of this 
anomaly, the transfer from ground to excited state of N 2 should have been changed by a factor of 34 also. Since all 
the spectra are essentially identical (aside from different levels of instrument noise) we determined that the 
population distribution is not related to the laser intensity. Thus, SRP is not significant and what we are observing is 
vibrational non-equilibrium unrelated to the CARS measurement process. 



Figure 3. CARS spectra acquired using green laser energies differing by a factor of 34, which was 
accomplished by rotating the angle of a half wave plate in combination with a polarizer used to attenuate the 
beam. 



2120 2170 2220 2270 2320 

Wave number, cm’^ 

Figure 4. CARS spectra acquired near the center of the duct at x/h = 55.6 (Configuration C), for constant 
facility stagnation temperature and pressure but varying mole fractions of steam. 

During the Configuration C test campaign (with essentially identical CARS instrumentation to A) we acquired 
CARS data sets at the entrance plane with various levels of steam introduced in the facility heater. Since steam is 
known to greatly increase the rate at which N 2 relaxes to thermal equilibrium,^^ this experiment was expected to 
confirm the hypothesis that non-equilibrium was a result of processes in the flow and provide useful information on 
relaxation rates. Figure 4 shows experimental CARS spectra (200 shot averages) acquired near the center of the duct 
for cases with no steam (i.e., “clean air”), and 1.9%, 3.7% and 5.7% steam. The plots are normalized to the peak 
signal to overlap the rotational band structure of the ground state of N 2 for each spectrum. The rotational band 
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structure of the ground state of the N 2 and the O 2 are seen to be identical in all cases, but the height of the first 
excited vibrational band is reduced by steam, as expected. The spectra for 3.7% and 5.7% steam are almost the 
same, and we assume they are at the equilibrium condition. For 1.9% the excited band of N 2 is a little high, but for 
O 2 it is equilibrium. 




Figure 5. CARS measured spectra and theoretical fits for clean air (left) and 5.7% steam (right) at the center 
of the plane for Configuration C. The theory for clean air assumed Ty^m = 1200 K, Tv,o 2 = 1000 K, and Tj. = 735 
K. The theory for 5.7% steam assumed Ty=Ty = 735 K. 

Figure 5 shows fits to theory for these same spectra with no steam and 5.7% steam using CARSFT. The fits were 
performed by trial and error variation of rotational and vibrational temperatures, using the known pressure and 
composition. (This fit with no steam is comparable in appearance and gives the same result within uncertainty to the 
computed fit described in the next section.) In the case of 5.7% steam, the best fit occurred with equal rotational and 
vibrational temperatures {Ty = 735 K ± 25 K), i.e., with equilibrium, whereas for no steam, different vibrational 
temperatures for N 2 and O 2 were found {Tr = 735 ± 25 K, = 1200 ± 40 K, Ty^o 2 = 1000 ± 40 K). Thus, the 
vibrational temperature of N 2 is frozen at or close to the temperature of the flow exiting the heater (prior to 
acceleration in the nozzle) while for O 2 it is frozen at a somewhat lower temperature. The uncertainties are 
estimates based again on trial and error and on the sensitivity of the fit to the temperature parameters. 

B. Spatial Variation of and Ty for Clean Air 

The measurement planes for both configuration A and C were mapped by acquiring 200 spectra at each of about 
170 points on a regular 3 mm (in z) x 2.4 mm (iny) spaced rectangular grid. Some of these grid points lay within 
the window slots, which were ventilated by cooling air and so registered lower temperatures. In this analysis, the 
spectra were averaged before fitting. Our library fitting routine is not adapted to fitting vibrational temperature 
separately from rotational temperature. Our strategy was thus to fit the complete data sets several times with 
different fixed vibrational temperatures and compare the residuals for each fit. For Configuration C, six libraries 
were generated, with Ty^m = Ty 02 of 500, 700, 900, 1100, 1300 and 1500 K. To find the vibrational temperature of 
N 2 at each point, the data set was fitted with each library, blocking from the fits the portion of the spectrum 
containing the O 2 Q-branch (from 2120 cm'^ to 2196 cm'^). The fit residual at each point in the grid was fit as a 
function of temperature to a cubic spline curve and the minimum of this curve was the experimental Ty^m- To find 
the vibrational temperature of O 2 at each point, the procedure was repeated except that the portion of the spectrum 
containing the N 2 Q-branch was blocked (2265 cm'^ to 2352 cm'^). The best rotational temperature at each point 
was found by generating an additional library using the average fitted (1200 K) and Ty^o 2 (1000 K), and fitting 
the spectrum for T^. This was done three times; first using the complete spectrum, second blocking the N 2 Q-branch, 
and third blocking the O 2 Q-branch. When the N 2 was blocked the temperature was 21 K higher, on average, 
whereas when the O 2 was blocked the temperature was 1 1 K lower. The results with no blocking are presented and 
the spread in temperature is taken as an indication of uncertainty. A similar approach was adopted in analyzing the 
Configuration A data, differing in some minor details because of differences in the experiment. 

The spatial distributions of Ty^N2-> Tv, 02 and Ty are plotted in Fig. 6; the left and right hand column in the figure 
contain the Configuration A and C results respectively, while the first, second and third rows contain Ty^N 2 , Ty 02 and 
Tr respectively. The measurements of Ty^o 2 for Configuration A are not shown because the precision errors are too 
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high to yield meaningful results. All these plots show a deficit in temperature in the top right hand quadrant of the 
flow. This deficit was previously reported in Ref 14 for Configuration A, but it is significant that it is still present at 
the further downstream location in Configuration C. Since this deficit occurs even in the vibrational temperature of 
N 2 -which is frozen at close to the facility stagnation condition-the deficit must have existed in the flow before 
entering the nozzle and be due to non-uniformity in electrical heating of the air. The plot also shows bands of lower 
temperature at the side walls at z/h = ±3 caused by the cooling air entering from the window slot. As will be 
discussed in the next section, the side wall boundary layers are visible in the range 2 < \z/h\ <2.75, manifest by an 
increased rotational temperature (due to viscous heating) and slightly decreased vibrational temperature of N 2 . The 
top and bottom wall boundary layer cannot be seen due to lack of data near enough to the walls, which are at z/h = 
±2. 




z/h z/h 

Figure 6 . Maps of temperature measured by CARS. The figures on the left are Configuration A, on the right 
Configuration C. From top to bottom, Tv,n29 ^v,o 2? and The duct side walls are at z/h = ± 3 . Color scales 
indicated in the legends on the right are consistent left to right but different top to bottom. 
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C. Comparisons of Measured Temperatures with Calculations 

In this section, the results of the computational fluid dynamics calculations using the LAURA code are compared 
to the CARS measurements. The non-equilibrium simulation computes a translational temperature (equal to the 
rotational temperature) and a single vibrational temperature, T^. To enable the comparisons, a single lumped 
temperature T^jump was calculated from the CARS individually measured N 2 and O 2 vibrational temperatures by 
solving Eq. (1). 


0.79 Ev(0y^^2,Tyjump) + 0-21 Ev{6y^o2,Ty iufnp) ~ 0-79 Ev(0yj^2,Eym) + 0.21 Ev(0y O2,Ty O2) (1) 

R6 

We define Ev = , where R is the universal gas constant, and 6y is the characteristic vibrational temperature 

exp (f) 

(3390 K for N 2 and 2270 K for 62 ).^’^^ For example, given Ty^^2 = 1200 K, and Ty 02 = 1000 K, the solution is Tyj^mp 
= 1148K. 

The temperatures calculated by the LAURA code, assuming non-equilibrium, are plotted against y/h in Fig. 7, 
with (a) corresponding to the Configuration A CARS measurement plane and (b) to the Configuration C CARS 
plane. The lumped vibrational temperature in the freestream is 1140 K for A and 1115 K for C, while the 
rotational/translational temperature is 662 K for A and 691 K for C. The decrease in vibrational temperature is due 
to vibrational relaxation while the increase in rotational temperature is mostly due to the growth of boundary layers 
on the top and bottom wall of the isolator duct, which reduces the Mach number in the freestream, raising 
temperature and pressure. The calculation is likely to underestimate the actual compression by the boundary layers 
since it is a 2D calculation while the flow path is rectangular in cross-section (aspect ratio of 1.5) with side wall 
boundary layers not captured by the calculation. The top and bottom wall boundary layers are visible in the 
calculation as an increase in rotational temperature approaching the walls and a decrease in vibrational temperature. 
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Figure 7. Profiles of temperature on the center plane (z = 0). Comparisons between measured vibrational and 
rotational temperatures and computational result using the LAURA code, (a) Configuration A (x/h = 13.6), 
(b) Configuration C (x/h = 55.6). 


The CARS measured temperatures along the grid line closest to z = 0 are also shown in Fig. 7, for comparison 
with the calculation. For Configuration A, despite the previously mentioned spatial non-uniformity in temperature 
from the heater, the agreement of the measured rotational temperature with calculation is quite good, while the 
measured TvN 2 is greater than the calculated Tyj^mp, as expected (since Ty 02 is less). Notice that Tv,n 2 > 1200 K where 
y < 0, meaning that, in this part of the flow, the temperature entering the nozzle is higher than the nominal facility 
heater temperature (which was 1200K). For Configuration C, the measured rotational and lumped vibrational 
temperatures are, on average, a little higher than calculated. This difference is due to the spatial non-uniformity in 
heater temperature and the growth of the side wall boundary layers, one effect of which is to increase the overall 
pressure and temperature above that of a 2D calculation. Like the calculation, the CARS vibrational temperature of 
N 2 is a little less at C than at A due to vibrational relaxation. The calculation shows an increase in Ttrans above the 
freestream value in the top and bottom wall boundary layers, but a decrease in Tyjump- As discussed in the previous 
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section, this behavior is not seen in the CARS data (shown in Fig. 6) near the top and bottom walls due to lack of 
measurements there, but something similar is seen near the side walls. 

D. Effect of Steam on Wall Pressure and Flow Temperature 

Figure 8 shows measured and calculated wall static pressures, normalized by the facility heater pressure and 
plotted as a function of distance down the duct {x/h). The plot includes measurements with clean air and steam, as 
well as calculations assuming non-equilibrium and equilibrium. As the distance increases, the calculated pressure 
drops within the nozzle to a minimum at the nozzle exit, then rises gradually in the isolator due to growth of top and 
bottom wall boundary layers. For clean air, the measured pressure at the first point (after the nozzle exit) is in 
agreement with the non-equilibrium calculation, as is the general trend of increasing pressure in the isolator. 
However, there are oscillations in the measurements not observed in the calculations, presumably due to reflecting 
shock waves. The average slope of the measurement is a little higher than the slope of the calculation due to side 
wall boundary layers not reproduced by the 2D calculation. The effect of steam is to raise the pressure at the nozzle 
exit by from 4.5% (for 1.9% steam) to 6.4% (for 5.7% steam) and this difference is maintained throughout the 
isolator. Part of this change is due to a change in thermodynamic properties, but this effect is small as evidenced by 
the small change between the measurements as the percentage of steam is varied. On the other hand, the change due 
to steam is reproduced by switching the calculation from non-equilibrium to thermal equilibrium, consistent with the 
CARS observation that addition of steam brings the system towards equilibrium. 



Figure 8. Pressure on the interior flow surface of the flow path (z = 0): measured pressure for clean air and 
several steam levels, and computational results for non-equilibrium and equilibrium. 
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Figure 9 shows a map of the CARS temperature for Configuration C, 1.9% steam, spectra fitted assuming rv,N 2 = 
950 K and thermal equilibrium of O 2 . The temperature distribution in the map is generally similar to the distribution 
of rotational temperature for Configuration C, clean air shown in Fig. 6, although temperatures are on average about 
44 K (about 6%) higher. The CARS results at z = 0 are also shown in Fig. 10 as a function of y/h and compared to 
the measurements for clean air and the calculations. The measured increase in temperature from clean air to 1 .9% 
steam is accurately reproduced by the calculated increase from non-equilibrium to equilibrium, although the overall 
level of the calculation is low due to the neglect of the side wall boundary layers. 



Figure 10. Profiles of temperature on the center plane for Configuration C. CARS measurements with clean 
air and 1.9% steam. Computational results for non-equilibrium and equilibrium. 

Discussion and Conclusions 

In this work, we have measured the vibrational and rotational temperatures in two planes in a constant area duct 
mounted to the exit of a Mach 2 nozzle at the University of Virginia Supersonic Combustion Facility. With clean air 
flowing in the facility at a nominal stagnation temperature of 1200 K, we have measured non-equilibrium vibrational 
temperatures of N 2 and O 2 . Indeed, the N 2 vibrational temperature of about 1200 K was close to the facility heater 
stagnation temperature, while the vibrational temperature of O 2 was little lower (1000 K). This result was at first a 
surprise since, although vibrational freezing is commonly encountered in hypersonic nozzles, it is not typically 
expected in supersonic ones. The existence of vibrational freezing is a function of the particular operating 
conditions. Supersonic wind tunnels are not usually operated with heated air while hypersonic wind tunnels 
typically are, so vibrational modes of energy are simply not populated for supersonic nozzle flows. On the other 
hand, since the time required to establish thermal equilibrium grows steeply with decreasing temperature (and 
pressure),^ for a higher heater temperature, the air might not drop below the threshold in temperature at which 
freezing occurs. The relatively small scale of the facility is also a factor. Computational fluid dynamics calculations 
using established models for vibrational relaxation were able to predict the rotational and “lumped” vibrational 
temperatures. Measurements were also made in the facility with several molar percentages of steam (1.9%, 3.7% 
and 5.7%). Steam was found to bring the flow to thermal equilibrium, consistent with expectations based on the 
effect of steam on relaxation rates known from the literature. Steam significantly increases both temperature and 
pressure at the nozzle exit due to the change from vibrationally frozen to thermal equilibrium flow. This work also 
confirms previous observations of non-uniformities of temperature entering the facility nozzle, and shows that this 
non-uniformity persists at least to the end of the facility isolator. Both vibrational non-equilibrium and temperature 
non-uniformity are expected to affect the distribution of pressure and combustion in the scramjet flow path. This 
greater understanding of the flow in the nozzle and isolator will substantially assist our interpretation of 
measurements acquired in the dual-mode scramjet combustor experiments being conducted in the facility. 
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